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Abstract 

An analysis of traveling wave solutions of pure cross-diffusion systems, i.e., sys- 
tems that lack reaction and self-diffusion terms, is presented. Using the qualitative 
theory of phase plane analysis the conditions for existence of different types of wave 
solutions are formulated. In particular, it is shown that family of wave trains is a 
generic phenomenon in pure cross-diffusion systems. The results can be used for 
i-^p ■ construction and analysis of different mathematical models describing systems with 

O^. directional movement. 
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^ : 1 Introduction 

O ■ Using so-called cross-diffusion terms in mathematical models in natural sciences, in partic- 

ular in biology, is nowadays a routine trick to incorporate into classical models additional 

K^ ! effects such as taxis movement, crowding, pursuit, or evasion (for a review see [18f|). Two 

^ [ particular examples, which received especial attention amongst other models, are the 

Keller-Segel model of chemotactic movement [lO|, 112, ll7| {taxis is defined as behavioral 
response to a directional stimulus), and the Lotka-Volterra type cross- diffusion model to 
describe predator pursuit of an evading prey [7], Isj, [uJ, [l3| . 
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A distinct feature of many studies on cross-diffusion models (PDE systems where 
diffusion matrix is not diagonal) is that usually a particular model is analyzed, where terms 
of reaction and diffusion coefficients are given in fixed explicit form; notable exceptions 
are given by |3|, |5|, lid]. It also appears that almost all the models considered contain 
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reaction terms, only one exception we know of notwithstanding [8j, although it might be 
argued that systems without reaction terms possess numerous solutions which may be 
relevant in biological modeling. 

Taking a step further, in this note we consider the systems that lack the self- diffusion 
terms; PDE systems we analyze possess only cross-diffusion (as, e.g., in p, [iJ], but with 
nonlinear cross-diffusion coefficients). We dub such systems without reaction and self- 
diffusion terms pure cross- diffusion systems. In particular, we study the PDE system in 
the following form: 

ut = di{f{u,v)v^)^, 
vt = d2{g{u,v)ux)x, 

where u = u{x,t), v = v{x,t), — oo < x < oo, t ^ 0, di,d2 are real constants, and 
f{u,v), g{u,v) are arbitrary functions whose properties will be specified later. Qualita- 
tively system ([T]) describes mutually dependent movement of conserved entities u and v 
on gradients of their respective counterparts. Of particular interest for us are the travel- 
ing wave solutions of system (Q, which, e.g., can describe a replacement process of one 
species by another species. As usual, a wave solution of ([T]) is a bounded solution having 
the form u{x, t) = u{x + ct) = u{^), v{x, t) = v{x + ct) = v{C,), where c is the wave speed. 
On substituting these traveling wave forms into ([T]) we obtain the ODE wave system 

cu^ = di{f{u,v)v^)^, 
cv^ = d2{g{u,v)u^)^. 

Note that we can integrate system ([2]), and therefore, (E]) is essentially two dimensional, as 
the original PDE system. Each wave solution of ([T]) has its counterpart as a bounded orbit 
of ([2]), which, due to the particular structure of ([2]), is an orbit on a phase plane. There is 
a known correspondence between wave solutions of ([T]) and bounded orbits of ([2]): ^ wave 
front in, e.g., u component corresponds to a heteroclinic orbit of ([2]) that connects singular 
points with different u coordinates; A wave impulse in, e.g., u component corresponds to 
a heteroclinic orbit of ([I]) that connects singular points with identical u coordinates or 
to a homoclinic curve of a singular point; A wave train corresponds to a closed orbit of 
system ((21). 

In the next section we show that, under some natural restrictions on functions f{u,v) 
and g{u,v), all possible traveling wave solutions of ([T]) can be classified. In particular, a 
family of wave trains of ([T]) is a generic phenomenon intrinsic to ([T]). Section 3 provides 
some examples and extensions. Section 4 is devoted to discussion and conclusions. 

2 Wave solutions of the separable models 

In this note we confine ourselves to a particular case of ([I]), whereby the functions f{u, v) 
and g{u,v) can be represented as products of functions that depend on one variable: 

f{u,v) = fi{u)f2{v), g{u,v) = gi{u)g2{v). (3) 



Here the functions fi{x), gi{x), i = 1,2 are smooth or rational. The models of the form 
([T]) for which (jS]) holds will be termed separable pure cross diffusion models, or separable 
models for short. 

The rationale behind the choice of the from of the functions / and (7 as in ([3]) is 
twofold. Firstly, such assumption allows complete analysis of possible structures of the 
corresponding wave system, and thus, exhaustive classification of possible wave solutions 
of ([T]) as was done for different models in [4, |5|. Secondly, such forms are abundant in the 
modeling literature being the consequence of a biased random walk on the lattice (see, 
e.g., 0). 

For the separable model and after integration and rearrangement wave system ([2]) 
takes the form 

u^ = cd2{v-ri)/g{u,v), 
v^ = cdi{u- fi)/f{u,v), 

where cq, cfi are constants of integration that in general depend on boundarv conditions, 
and which are considered here as new parameters of the system (see also |5|). After the 
change of independent variable 



did2f{u,v)g{u,v) 



which is well defined and smooth for any {u,v) except where f{u,v)g{u,v) = or where 
f{u,v) or g{u,v) are not defined (see [ij), we obtain 

H = difi{u)f2{v){v-ri), 
V( = d2gi{u)g2{v){u- ^). 

In the case if f{u,v) or g{u,v) are rational, the change of independent variable is chosen 
such that the resulting wave system is well defined for all u, v (for examples see the next 
section). Therefore, without loss of generality, we assume that our wave system has the 
form 

u^ = di<^i{u)<^2{v){v-ri), 
V( = d2i'i{u)ijj2{v){u - n). 

which is defined for all m, v. 

First we note that system (Q is integrable. Function 

(v — rj) — -— dv — d2 (u — /i) — -— du + const 
i^2{v) J (pi{u) 

represents its integral. We remark that under the assumption (pi{u) = const, ^2(^) = 
const system ([6]) is hamiltonian, and H{u, v) is its Hamiltonian, which automatically 
yields that the singular points of (E]) can be only topological saddles or centers. 



We assume that the following conditions of non- degeneracy are fulfilled (these condi- 
tions are natural as far as we consider 77, /i, c the only parameters in the system ([H])): 

(CI) V^i(m) and ipi{u) have no common roots; 
(C2) f2{v) and il'2{v) have no common roots; 
(C3) fi{x) and ipi{x), i = 1,2, have no multiple roots. 

System ([6]) always has singular point (/i, 77) . Coordinates of other singular points are 
{n,Vi), {iii, rj), {ui,Vi), where Ui and Vi are the solutions of one of the systems: 

ipi{u)=Q, Mv) = Q, (7) 



or 



^2{v) = 0, Mu) = 0. (8) 

To infer the possible types of the singular points of ([U]) consider the Jacobi matrix of 

m-. DF{u,v) = {a,,)2^2.tJ = l,2: 

an = di{v - ri)^2{v){<^i{u))'^, 

a\2 = di!fi{u){{v - r])[}p2[v))'^ + (^2(f)), ,g. 

^21 = d24)2{v){{u - /i)(^i(M))^ + i)\{u)), 

a-22 = d2{u - ^i)il)i{u){il)2{v))'^. 

Direct substitution of singular point coordinates into ([9]) and evaluation of trace and 
determinant of DF implies the following proposition. 

Proposition 1. Let us assume that conditions (C1)-(C3) hold. We also assume that Ui 
{vi) do not coincide with fi {rj, respectively). Then 

1. Singular points of ([6]) satisfying (JTj) or ([H]) can be saddles or nodes only; 

2. Singular points of ^ where one or both coordinates are rj or fi can be centers or 
saddles only. 

Remark 1. If we relax condition (C3) retaining (CI) and (C2) then we have to add to 
the list of possible singular points of ([U]) saddle-node points. 

Remark 2. To prove that singular points that have eigenvalues with zero real part 
in linear approximation are indeed centers we can use integrability of system ([U]) and 
Theorem 2 in [2, p. 75]. 

From Proposition [1], integrability of system ([6]), and special structure of the wave 
system ([6]) follows that the structure of the phase plane, and, consequently, of the bounded 
orbits of (Q are completely determined by the types and mutual positions of singular 
points of ([H]). The phase plane is divided into orbit cells that represent either a rectangular 
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Figure 1: A typical structure of the phase plane of the wave system of the separable 
model ([T]). Here the bold dots represent singular points of i^, the arrows show direction 
of the vector field. The plane is divided into orbit cells that contain families of bounded 
orbits: / is a cell that corresponds to impulse-front solutions; // is a cell that corresponds 
to front-front solutions; /// is a cell that corresponds to traveling train solutions, the 
homoclinic surrounding the closed orbits corresponds to impulse-impulse solutions 



cell with orbit structure determined by the types of the corner singular points (as in [5i] 
or a cell filled with closed orbits surrounding singular point of the center type (an example 
of such phase plane is given in Fig. [1]). The former case corresponds to fronts or impulses 
of the initial system ([T]) , and the latter to the wave train solutions of system ([T]) . 
Summarizing, we obtain the basic result of the present note: 

Theorem 1. The separable pure cross- diffusion system ([T]), satisfying (C1)-(C3) can 
only possess traveling wave solutions of the following kinds: i. impulse-front solutions; ii. 
front-front solutions; Hi. impulse-impulse solutions; iv. family of wave train solutions. 

The family of wave train solutions exits in domain Q if f{u, v) and g{u, v) are analytical 
in Q and did2f{u, v)g{u, f ) < for {u, v) G Q. 



3 Examples 

Example 1. The simplest possible separable model ([T]) is obtained if we put 

fiiu) = f,{v) = g,iu) = g2iv) = I. (10) 



The integral in this case is given by 

c?2(m — /i)^ — di{v — rff = const. 

From which it follows that system ([T]) with ( TTOl) has no traveling wave solution if did2 > 
and has a family of traveling trains ii did2 < 0. 

Remark. Note that in case did2 < system ([T]) is equivalent to the beam equation 
Utt = —Duxxxx, D > 0, which, how it is well known, possesses wave solutions. 



Example 2. Let 



f2{v) = gi{u) = const 

h{u) = {u-ar,g2{v) = {v-br,n,meK ^ ' 

It is straightforward to prove 

Proposition 2. Separable model ([T]) with fill I) has a family of wave trains if did2 < and 
has a front-front solution, that corresponds to a heteroclinic connecting node and saddle, 
if did2 > 0. 

Example 3. Let 

f,{u)=u\f2{v) = ^, 

1 ' (12) 

giiu) = —, g2iv) =v"', 

where m,n G N. After a suitably chosen change of independent variable in the wave 
system (^ we obtain 



U( = (iiM^"(f — 7]), 

f^ = d2v'^"^{u - jj). 



2^:^. .. (13) 



Having in mind biological applications of the pure cross-diffusion models we assume that 
yU > 0, r] > and analyze the phase plane in the first quadrant. System ( TT3l) has two 
singular points: (0,0) and (/i, 77). The point (/i, r]) can be a saddle (if d\d2 > 0) or a 
center {d\d2 < 0). If we restrict our attention only to the cases /i > 0, r^ > 0, m, f > 0, 
then traveling wave solutions of ([T]) with (IT3l) are exactly the same as in Proposition [2l 

Let us illustrate this case. In Fig. [2] the first quadrant of the phase plane of the 
wave system of (fT]).( IT2l) is shown for two topologically nonequivalent cases. The first case 
corresponds to the wave train in the original PDE system, and the second case corresponds 
to front-front solution. In Fig. [3] numerically obtained solutions of flTl). flT^ are presented. 
To approximate cross-diffusion term an "upwind" explicit scheme was used J15l| which is 



often used for cross-diffusion systems (e.g., |l9f|). 





Figure 2: The structure of the first quadrant of the phase plane of the wave system ([2]), 
(fT2l) . In case did2 < a family of closed orbits can be found, which corresponds to the 
wave trains; in case did2 > there is only one bounded trajectory that connect the origin 
and the saddle point {fi, rj) 

Example 4. We can use obtained results to prove existence of solutions of a special 
form when simple reaction terms are presented. For example, let us consider the following 
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Figure 3: Numerically obtained wave train solutions of system (ITl).( IT2l) at time moments 
to = O(bold curves) < ti < t2 < ts through equal time intervals. The arrow shows the 
direction of the waves 



system: 

ut = biu + diiu^/v'^v^)^, 
vt = b2V + d2{v"^/u''u^)^. 

If the constant bi, 62, n, m satisfy the following equality: 

61(1 -n) = 62(1 -m), 

then by the change 

u{x,t) = e''^^z{z,t), v{x,t) = e^^^w{x,t) 

system ^\M can be reduced to system ([T]), (fT2l) . Thus system ^\M has wave train solutions 
with amplitude changing with time. It is interesting to note that for some models of the 
form ^\M an exact solution can be found. 

Example 5. Consider the system 

Ut = U + diVxx, 



Vt = -IV + d2Uxx- 



Proposition 3. For 7 < 1 and di + d2 > system ( 1151) has a family of wave trains. 

Proof. We will prove proposition by presenting a solution of (1151) in explicit form. 

Let m(^) = cosfc^, f(^) = sin/c^, ^ = x + ct. Then system ( ITSl) transforms to the 
equations 

—ck sin k^ = sin k^ — dik'^ sin k^, 
ck cos k^ = —7 cos kC, — ^2^^ cos k^, 

from which 

dik^ -ck -1=0, 
d2k'^ + cA; + 7 = 

follows. For fixed di, d2, 7 the solution of the last system is given by 



_ yi -7 4 + di-f 



di + d2' v/(c/i + rf2)(l-7)' 

which is defined for 7 < 1 and di + d2> 0, which completes the proof. D 



Example 6. An obvious question is what kind of solutions can appear in pure cross- 
diffusion models if we do not demand that the functions f{u,v) and g{u,v) can be rep- 
resented as the product of functions of one variable. Despite the fact that the separable 
models possess very reach structure of possible traveling wave solutions, simple extension 
of the models considered can be shown to yield a different type of behavior. 
Consider the model ([T]) with the functions given by 

f{u,v) = — , 

, , V (16) 

9{u,v) = — — — , 

u + pv 

where a and /3 are positive constant. 

The wave system, after suitable change of independent variable, reads 

u^ = di{v-r])u{au + ri), 
V(^ = d2{u — ji)v{u + l3v). 

Assume that all the parameters, except for rfi, d2 are positive. Singular point (/i, rf) is 
a saddle, if did2 > and a center at linear approximation for did2 < 0. However, contrast 
to model in Example 3, and due to the fact that the wave system is now non-integrable, 
we cannot make conclusion on the type of the singular point in the case did2 < 0. If we 
calculate the first Lyapunov value h [2], we obtain h = did2^J-ri{a(3 — l){a^'^ — (3rj^), which 
shows that for parameter values where /i 7^ this singular point can be either stable or 
unstable nonlinear focus. 

Numerical experiments show that at least for some parameter values there is a limit 
cycle surrounding the singular point (see Fig. Hj). Thus, in this example, there is a wave 
train solution of the initial pure cross- diffusion system. 

4 Conclusions 

In the present note we have shown that pure cross-diffusion models possess very rich 
structure of possible traveling wave solutions even in a rather simplified, though still 
realistic, case of separable models, when the nonlinear cross-diffusion coefficients in system 
^ can be represented by product of functions depending on one variable. 

A distinct feature of pure cross diffusion models is the presence, under some additional 
condition, of the family of wave train solutions (such solutions were not found for several 
classed of cross-diffusion models describing taxis, see [J, Isj). As numerical experiments 
show, such wave trains can be observed at least for some time interval. Together with fam- 
ilies of wave trains, other typical cases of traveling waves can be observed. In particular, 
impulse- front, front-front, and impulse- impulse solutions are also a possibility. 

We emphasize that the separable model ([T]), even when the values of parameters /i, 77, c 
are fixed, possesses, in general, a family of traveling wave solutions; there are infinitely 
many bounded orbits of ([2]) that correspond to traveling wave solutions. 




Figure 4: Limit cycle in the wave system (TT7|) . Parameter values di = —1, d2 
0.7, /5 = 0.2, /i = 0.3, ?7 = 0.5 
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The presented results can be used for construction and analysis of different mathemat- 
ical models describing systems where directional movement is of particular importance. 
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